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dependent lysis, because the calcium chelator 
EGTA, which blocks the perforin pathway, also in- 
hibits up-regulation of FasL, as required in our stud- 
ies with primary human T cells. 
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addition of blocking antibody to TNF-a did not inhibit 
the cytotoxicity. 
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CD1+ macrophages, which were infected with live M. 
tuberculosis, were labeled with 100 wCi of ©'Cr (ICN, 
Costa Mesa, CA) for 1 hour and plated in a 96-well 
V-bottom plate at a final concentration of 4000 targets 
per 100 wl. Appropriate samples were incubated with 
blocking antibodies to CD1a, CD1b, or CD1c for 30 min 
before addition of the T cells. After a 9-hour incubation, 
target cell lysis was calculated by quantifying "Cr re- 
lease in a gamma-counter. The data are given as per- 
cent specific lysis, calculated as [(cpm release from ex- 
perimental — com spontaneous release)/(maximal re- 
lease — com spontaneous release) X 100]. The spon- 
taneous °'Cr release by macrophages in the absence of 
T cells was <15%. 

CD1* macrophages were pulsed with soluble M. 
tuberculosis extract (5 wg/ml) overnight, detached 
with 1 mM EDTA, and labeled with 100 Ci of ©'Cr 
for 1 hour. For inhibition of the interaction between 
FasL and Fas, the assay was done in the presence 
of blocking antibodies to FasL (27) (Pharmingen, 
San Diego, CA) or Fas (22) (Immunotech, West- 
brook). Degranulation of the cytotoxic granules 
was induced by initial treatment of T cells with 25 
mM Sr?* (Aldrich, Milwaukee, WI) for 18 hours. 
51Cr release was determined after a 4-hour incu- 
bation. Expression of Fas on the target cells was 
not affected by infection with M. tuberculosis, as 
determined by flow cytometry. 

CTLs (5 X 105) were incubated in the presence of 25 
mM Sr?* for 10 hours in a final volume of 1.5 ml. The 
amount of Na-CBZ-, -lysine thiobenzyl (BLT )-ester- 
ase in the supernatant was determined by the BLT- 
esterase assay (23). The supernatants (20 ul) were 
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of 1 mM 5,5’-dithio-bis-(2-nitrobenzoic acid) (Sig- 
ma), and 10 pl of a 0.1% Triton X-100 (Sigma) solu- 
tion. After a 30-min incubation at 37°C, the absor- 
bance at 405 nm was determined. 
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Mitochondrial DNA control region sequences were analyzed from 162 wolves at 27 
localities worldwide and from 140 domestic dogs representing 67 breeds. Sequences 
from both dogs and wolves showed considerable diversity and supported the hypothesis 
that wolves were the ancestors of dogs. Most dog sequences belonged to a divergent 
monophyletic clade sharing no sequences with wolves. The sequence divergence within 
this clade suggested that dogs originated more than 100,000 years before the present. 
Associations of dog haplotypes with other wolf lineages indicated episodes of admixture 
between wolves and dogs. Repeated genetic exchange between dog and wolf popu- 
lations may have been an important source of variation for artificial selection. 


The archaeological record cannot resolve 
whether domestic dogs originated from a 
single wolf population or arose from multi- 
ple populations at different times (1, 2). 
However, circumstantial evidence suggests 
that dogs may have diverse origins (3). Dur- 
ing most of the late Pleistocene, humans 
and wolves coexisted over a wide geograph- 
ic area (1), providing ample opportunity for 
independent domestication events and 
continued genetic exchange between 
wolves and dogs. The extreme phenotypic 
diversity of dogs, even during the early 
stages of domestication (1, 3, 4), also sug- 
gests a varied genetic heritage. Conse- 
quently, the genetic diversity of dogs may 
have been enriched by multiple founding 
events, possibly followed by occasional in- 
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terbreeding with wild wolf populations. 

We sequenced portions of the mito- 
chondrial DNA of wolves and domestic 
dogs. Initially, 261 base pairs (bp) of the 
left domain of the mitochondrial control 
region (5) were sequenced from 140 dogs 
representing 67 breeds and five cross- 
breeds and 162 wolves representing 27 
populations from throughout Europe, 
Asia, and North America (Fig. 1) (6). 
Because all wild species of the genus Canis 
can interbreed (7) and thus are potential 
ancestors of the domestic dog, five coyotes 
(Canis latrans) and two golden, two black- 
backed, and eight Simien jackals (C. au- 
reus, C. mesomelas, and C. simensis, re- 
spectively) were also sequenced. 

The control region of wolves and dogs 
was highly polymorphic (Fig. 1). We 
found 27 wolf haplotypes that differed on 
average by 5.31 + 0.11 (+SE) substitu- 
tions (2.10 + 0.04%), with a maximum of 
10 substitutions (3.95%). The distribution 
of wolf haplotypes demonstrated geo- 
graphic specificity, with most localities 
containing haplotypes unique to a partic- 
ular region (Fig. 1). Four haplotypes (W2, 
W7, W14, and W22) had a widespread 
distribution. In dogs, 26 haplotypes were 
found. Only haplotype D6 also occurred in 
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some gray wolves from western Russia and 
Romania (W6). Sequence divergence 
among dogs was similar to that found 
among wolves. Dog haplotypes differed by 
an average of 5.30 + 0.17 substitutions 
(2.06 + 0.07%), with a maximum diver- 
gence of 12 substitutions (4.67%). Mito- 
chondrial haplotype diversity in dogs 
could not be partitioned according to 
breeds. For example, in eight German 
shepherds examined, five distinct se- 
quences were found, and in six golden 
retrievers, four sequences were detected. 
Moreover, many breeds shared sequences 
with other breeds. For instance, dog hap- 
lotypes D4, D3, D5, and D1 were found in 
14, 14, 9, and 7 breeds, respectively. No 


Fig. 1. Substitutions and deletions (—) observed in 
261 bp of control region sequence from wolves 
W) and dogs (D). The dog sequence D13 had the 
same sequence as D4 except for an insertion of a 
67-bp tandem repeat. The numerals |, Il, Ill, and IV 
indicate assignments to the four clades of dog 
sequences. Wolf localities: Bulgaria (9 = 1, W7); 
Croatia (n = 5, W2); Estonia (n = 1, W10); France 
Nn = 2, W4); Finland (n = 2, W10); Greece (n = 7; 
W2, W5, W8, and W9); Italy (n = 12, W4); Poland 
n = 1, W8); Portugal (9 = 19; W1 and W2); 
Romania (n = 4; W5 and W6); Russia (n = 3; W6, 
W10, and W26); Spain (n = 46; W1 and W9); 
Sweden (n = 2; W2 and W10); Afghanistan (n = 3, 
W18); China (n = 3; W14, W19, and W27); India 
(n = 1, W12); Iran (n = 6; W16 and W17); Israel 
(n = 16, W11); Saudi Arabia (n = 7; W7, W12, 
W13, W14, and W15); Turkey (n = 2, W2); Alaska 
(n = 3, W20); Alberta (n = 1, W22); Labrador (n = 
8, W22); Mexico (n = 5, W25); Montana (n = 1, 
W22); Northwest Territories (9 = 3, W22); and 
Yukon (n = 3; W21, W238, and W24). Dog breeds: 
basenji (n = 1, D2); basset (n = 1, D6); boxer (n = 
1, D4); Norwegian buhund (n = 1, D1); bulldog 
n = 1, D6); Chinese crested (n = 2; D2 and D25); 
chow chow (n = 3; D1, D2, and D9); collie (n = 1, 
D1); border collie (9 = 3; D1 and D8); wirehaired 
dachshund (n = 3; D5 and D10); Australian dingo 
in = 4, D18); grey Norwegian elkhound (n = 9; D3 
and D8); Eskimo dog (n = 1, D23); German shep- 
herd (n = 8; D4, D5, D6, D7, and D19); greyhound 
n = 1, D9); groenendael (n = 1, D6); Mexican 
hairless (n = 6; D8, D6, D21, and D26); hamilton- 
stévare (n = 1, D5); Afghanistan hound (n = 3, 
D6); Alaskan husky (n = 2; D4 and D7); Siberian 
husky (n = 3; D8, D7, and D18); jamthund (n = 3; 
D7 and D8); keeshond (n = 1, D5); kuvasz (n = 


ane 


dog sequence differed from any wolf se- 
quence by more than 12. substitutions, 
whereas dogs differed from coyotes and 
jackals by at least 20 substitutions and two 
insertions. These results clearly support 
wolf ancestry for dogs. However, because 
mitochondrial DNA is maternally inherit- 
ed, interbreeding between female dogs and 
male coyotes or jackals would not be de- 
tected. More limited studies of nuclear 
markers support the conclusion that the 
wolf was the ancestor of the domestic dog 
(8). 

Several methods of phylogenetic anal- 
ysis, including maximum likelihood (9), 
maximum parsimony (10), minimum 
spanning networks (11), and_ statistical 
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D4); Leonberger (n = 2; D1 and D4); Norwegian lundehund (n = 1, D16); Mareema (n = 1, D=6); 
Pyrenean mastiff (9 = 1, D11); Newfoundland (n = 1, D4); otter hound (n = 1, D6); papillon (n = 2; D3 
and D4); poodle (n = 1, D8); toy poodle (n = 1, D6); pug (n = 1, D26); Chesapeake Bay retriever (n = 
1, D198); flat-coated retriever (9 = 3; D4 and D10); golden retriever (n = 6; D4, D6, D15, and D24); 
abrador retriever (9 = 6; D4 and D12); Rhodesian ridgeback (n = 1, D26); rottweiler (n = 2, D8); 
Samoyede (n = 3; D1, D4, and D4); St. Bernard (n = 1, D9); schipperke (n = 1, D4); giant schnauzer (n = 
8; D4 and D7); miniature schnauzer (n = 1, D9); English setter (n = 4; D3 and D5); Irish setter (n = 3; D1 
and D9); New Guinea singing dog (n = 2, D18); shar (n = 1, D26); Icelandic sheepdog (n = 1, D3); Old 
English sheepdog (n = 1, D5); shiba inu (n = 1, D20); Cavalier King Charles spaniel (n = 1, D17); Irish 


water spaniel (n = 1, D6); springer spaniel (9 = 1, D8); Tibetan spaniel (n = 1, D6); spitz (n = 1, D22); 
Japanese spitz (n = 1, D8); airedale terrier (n = 1, D7); border terrier (n = 2, D3); fox terrier (n = 2; D3 
and D14); Norfolk terrier (n = 2, D4); West Highland terrier (n = 2, D7); Tibetan terrier (n = 2; D2 and D9); 
wachtelhund (n = 1, D5); whippet (n = 1, D3); Irish wolfhound (n = 2, D11); and crossbreeds (n = 5; D1, 


D3, D4, DS, and D18). 
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parsimony (12), were used to investigate 
relationships among sequences. All analy- 
ses supported a grouping of dog haplotypes 
into four distinct clades, although the to- 
pology within and among clades differed 
among trees (13). As exemplified by the 
neighbor-joining analysis (Fig. 2A), three 
of the four monophyletic clades defined a 
larger clade containing all but three dog 
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Fig. 2. (A) Neighbor-joining tree of wolf and dog 
haplotypes (D13 excluded; see Fig. 1) based on 
261 bp of control region sequence (77). (B) 
Neighbor-joining tree of 8 wolf and 15 dog ge- 
notypes based on 1030 bp of control region 
sequence. The suffixes a, b, and c after the hap- 
lotype labels were used to distinguish identical 
261-bp sequences that have different 1030-bp 
sequences. Bootstrap support is indicated at 
nodes if found in more than 50% of 10,000 boot- 
strap trees. 
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haplotypes and a subset of wolf haplotypes 
(W4 and W5). Clade I included 19 of the 
26 dog haplotypes. This group contained 
representatives of many common breeds as 
well as ancient breeds such as the dingo, 
New Guinea singing dog, African basenji, 
and greyhound (14). Clade II included dog 
haplotype D8, from two Scandinavian 
breeds (elkhound and jamthund), and was 
closely related to two wolf haplotypes 
found in Italy, France, Romania, and 
Greece (W4 and W5). Clade III contained 
three dog haplotypes (D7, D19, and D21) 
found in a variety of breeds such as the 
German shepherd, Siberian husky, and 
Mexican hairless. Finally, clade IV con- 
tained three haplotypes (D6, D10, and 
D24) that were identical or very similar to 
a wolf haplotype (W6) found in Romania 
and western Russia, which suggests recent 
hybridization between dogs and wolves. 
Many breeds contained representatives of 
more than one dog haplotype grouping 
(Fig. 1). 

Because the overall bootstrap support 
for many of the internodes in Fig. 2A was 
low, 1030 bp of the control region were 
sequenced for 24 canids, including repre- 
sentatives of the four dog clades (15). 
Although the association of clades was 
different, the analyses of the longer se- 
quences provided stronger support for the 
four monophyletic groupings of dog hap- 
lotypes (Fig. 2B) (13). A Wilcoxon 
signed-rank test was used to assess the 
monophyly of dog clades (16). Monophyly 
of all dog haplotypes can be rejected, and 
monophyly of clades I, II, and HI is mar- 
ginally rejected (P = 0.0004 and P = 
0.053, respectively). In both trees, dog 
haplotype clades I and IV are most closely 
related to wolf sequences from eastern Eu- 
rope (Greece, Italy, Romania, and western 
Russia). 

The coyote and wolf have a sequence 
divergence of 0.075 + 0.002 (17) and 
diverged about one million years ago, as 
estimated from the fossil record (18). Con- 
sequently, because the sequence diver- 
gence between the most different geno- 
types in clade I (the most diverse group of 
dog sequences) is no more than 0.010, this 
implies that dogs could have originated as 
much as 135,000 years ago (19). Although 
such estimates may be inflated by unob- 
served multiple substitutions at hypervari- 
able sites (20), the sequence divergence 
within clade I clearly implies an origin 
more ancient than the 14,000 years before 
the present suggested by the archaeologi- 
cal record (21). Nevertheless, bones of 
wolves have been found in association 
with those of hominids from as early as the 
middle Pleistocene, up to 400,000 years 
ago (1, 22). The ancient dates for domes- 
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tication based on the control region se- 
quences cannot be explained by the reten- 
tion of ancestral wolf lineages, because 
clade I is exclusively monophyletic with 
respect to dog sequences and thus the 
separation between dogs and wolves has 
been long enough for coalescence to have 
occurred. To explain the discrepancy in 
dates, we hypothesize that early domestic 
dogs may not have been morphologically 
distinct from their wild relatives. Con- 
ceivably, the change around 10,000 to 
15,000 years ago from nomadic hunter- 
gatherer societies to more sedentary agri- 
cultural population centers may have im- 
posed new selective regimes on dogs that 
resulted in marked phenotypic divergence 
from wild wolves (23). 

Although individual breeds show uni- 
formity with respect to behavior and mor- 
phology, most breeds show evidence of a 
genetically diverse heritage because they 
contain different haplotypes. Moreover, 
dog sequences cluster with different groups 
of wolf haplotypes. Therefore, after the 
origin of dogs from a wolf ancestor, dogs 
and wolves may have continued to ex- 
change genes. Backcrossing events could 
have provided part of the raw material for 
artificial selection and for the extraordi- 
nary degree of phenotypic diversity in the 
domestic dog. Domestic species of plants 
and animals whose wild progenitors are 
extinct cannot be enriched through peri- 
odic interbreeding, and change under ar- 
tificial selection may be more limited. 
Consequently, the preservation of wild 
progenitors may be a critical issue in the 
continued evolution of domestic plants 
and animals. 
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